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ABSTRACT 



During the formation of the large scale structure of the Universe, matter ac- 
cretes onto high density peaks. Accreting collisionless dark matter (DM) forms caus- 
ed ' tics around them, while accreting collisional baryonic matter (BM) forms accretion 

shocks. The properties of the accreting matter depend upon the power spectrum of 
the initial perturbations on a given scale as well as the background expansion in a given 
cosmological model. In this paper, we have calculated the accretion of DM particles in 
one-dimensional spherical geometry under various cosmological models including the 
Einstein-de Sitter universe, the open universe with Cl < 1, and the flat universe with 
Q\ = 1 — o . A density parameter in the range 0.1 < tt < 1 has been considered. 
. The initial perturbation characterized by a point mass at the origin has been consid- 

' ered. Since the accretion shock of BM is expected to form close to the first caustic of 

\Q \ DM, the properties of the accreting BM are common with those of the DM. Hence, 

ON ■ the accretion calculations with DM particles have been used to find the position and 

velocity of the accretion shock and the cluster mass inside it. The average temperature 
of BM has been estimated by adopting simplifying assumptions. The velocity of the 
accreting BM around clusters of a given temperature is smaller in a universe with 
smaller i7 Q , but only by up to ~ 24% in the models with 0.1 < G < 1. Thus, it would 
be difficult to use that quantity to discriminate among the cosmological models. How- 
ever, the accretion velocity around clusters of a given mass or a given radius depends 
more sensitively on the cosmological models. It is smaller in a universe with smaller fi 
by up to ~ 41% and ~ 65%, respectively. So, it can provide a better signature of the 
background expansion for different cosmological models. Although the existence of the 
fyN , caustics and the accretion shocks may not be confirmed by direct x-ray observations, 

■ the infalling warm gas of 10 4 — 10 5 K upstream of the shocks may be observed as the 

absorption systems of quasar emission lines. According to this study, the suggestion 
made by Kang, Ryu, & Jones (1996) that the large scale accretion shocks around 
clusters of galaxies can serve as possible acceleration sites of ultra high energy cosmic 
rays above 10 18 eV remains plausible in all viable cosmological models. 
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1 INTRODUCTION 

Recent satellite x-ray observations have made clusters of 
galaxies increasingly important to cosmology as a probe 
into the large scale structure of the universe. Being mas- 
sive and rare, their abundance in the local and distant uni- 
verse carries vital information on the initial density fluc- 
tuations and the matter content of the universe (Lubin et 
al. 1995; Eke, Cole, & Frenk 1996). Also being relatively 
young dynamically, the details of their structures can pro- 
vide us with some signatures left over from the formation 



epoch as well as information on the background cosmology 
(Crone, Evrard, & Richstone 1994; Navarro, Frenk, & White 
1995; Tsai & Buote 1995). A general consensus seems to be 
that the standard cold dark matter model (i.e. , Q = 1 and 
h — 0.5) normalized to the COBE DMR measurement of 
the anisotropies in the cosmic background radiation (i.e. , 
<T8 > 1) has serious difficulties in explaining the observed 
properties of x-ray clusters, such as the cluster abundance 
(Kang et al. 1994), the baryon fraction in clusters (Lubin et 
al. 1995), and the contribution of cluster emission to the x- 
ray background (Kang et al. 1994; Kitayama & Suto 1996). 
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On the other hand, an open CDM model or a flat CDM 
model with a cosmological constant with a smaller value 
of erg seems more consistent with many observations, and so 
has become popular recently (Cen & Ostriker 1994; Ostriker 
& Steinhardt 1995). 

Although the application of the Press-Schechter formal- 
ism (Press & Schechter 1974) to the hierarchical cluster- 
ing model of cluster formation has been successful in get- 
ting good agreement with the results of TV-body simulations 
(Lacey & Cole 1994), this semi-analytic approach cannot 
include the effects of non-equilibrium hydrodynamics that 
could be important during some phases in cluster evolution 
(Kang et al. 1994). Some studies suggest that the amplitude 
of the density power spectrum, its, can be constrained by the 
local cluster abundance (Eke, Cole, & Frenk 1996; Viana & 
Liddle 1995) and the range of the allowed value of Q can be 
narrowed down by examining the evolution of cluster abun- 
dance at low redshifts (i.e. , z < 1) (Bahcall & Cen 1992; 
White, Efstathiou, & Frenk 1993). Despite this, the obser- 
vational and theoretical/numerical errors involved in such 
procedures seem to be too big to make any consistent pre- 
dictions on those key cosmological parameters (Castander 
et al. 1995; Luppino & Gioia 1995; Tsai & Buote 1995). It is 
also unlikely that all crucial physics can be included in any 
variants of the treatments. 

The properties of x-ray clusters, other than the abun- 
dance and the baryon fraction mentioned above, which 
have been investigated for various cosmological models, in- 
clude the cluster-cluster correlation function (Bahcall & Cen 
1992), the density and velocity profiles (Crone, Evrard, & 
Richstone 1994), and the degree of substructures in the in- 
tracluster medium (ICM) (Tsai & Buote 1995). The lat- 
ter two, associated with the internal structure of clusters, 
are closely related with recent mergers and accretion onto 
the cluster mass scale (~ 8/i _1 Mpc), so dependences on 
Q, n (the power-law index of the initial power spectrum), 
and as could be degenerate. In addition to that, statisti- 
cal treatments of observed data of those quantities have a 
harder time discriminating clearly among different cosmo- 
logical models compared to the statistics of the cluster abun- 
dance. Even though there are tantalizing possibilities that 
the various properties of x-ray clusters can indeed be used 
to unveil the fundamental nature of the universe, further 
improvements by a factor of at least a few in both theoreti- 
cal and observational fronts seem to be required in order to 
make some solid predictions. 

In this paper, we have examined one more physical 
property associated with x-ray clusters, the accretion flow 
infalling toward the clusters. Matter accretes onto the high 
density peaks continuously throughout the history of the 
universe. Its rate on the scale of the cluster mass will be 
determined by the initial density fluctuations and the back- 
ground cosmology. Fluctuations continue to grow in the 
f2 = 1 universe, while they stop growing at a redshift 
z ~ 1 — 1 in a universe with small Q (Peebles 1993). On 
the other hand, the accretion on larger scales gets stronger 
for a smaller n for a scale-free power spectrum. Crone, 
Evrard, & Richstone (1994) examined the radial velocity 
profile of the accretion flow through JV-body simulations of 
different cosmological models with scale-free power spectra. 
They found that the density profile is flatter and the accre- 
tion regime is stronger for higher fi and smaller n. 



If the initial density perturbation is scale-free, the accre- 
tion of both baryonic matter and dark matter in the fi = 1 
universe can be described semi-analytically with a self- 
similar solution (Fillmore & Goldreich 1984; Bertschinger 
1985). But in a universe with fl D ^ 1, the self-similar solu- 
tion is not possible, since, in addition to the time scale mark- 
ing the transition of the initial perturbation to the nonlinear 
regime, another time scale such as the cosmic time t 1 / 2 when 
|1 — n| = 1/2 enters the problem. 

Here, we use a one-dimensional spherical iV-body code, 
which can follow the evolution of collisionless dark matter 
particles, to find the properties of accretion flows onto ob- 
jects of cluster mass scales in different cosmological mod- 
els: the Einstein-de Sitter universe with fi = 1, the low- 
density, open universe with Q a < 1, and the low-density, 
flat universe with Q < 1 and Qa = 1 — fi (from non-zero 
cosmological constant, A). While most previous studies on 
accretion focused on the density profile inside the collapsed 
objects developed from infall, our primary interests lie in the 
properties of spherical accretion flows outside the objects. In 
the accretion of both dark matter (DM) and baryonic mat- 
ter (BM), the collisionless DM forms caustics around the 
overdensity, while the collisional BM forms accretion shocks 
stopping the infalling material. In this situation, as pointed 
by Bertschinger (1985), the position of the accretion shock 
of the gas with 7 = 5/3 is very close to the position of the 
first caustic of the dark matter particles. In addition, since 
in the region outside the shock the accretion solution for BM 
is almost identical to that for DM, the infall velocity of BM 
upstream of the accretion shock can be described by that 
of the DM. Hence, hydrodynamic calculations to follow BM 
are not necessary if we are interested only in the position 
of the accretion shocks and the properties of the accretion 
flows. 

Of course, one-dimensional treatments cannot include 
important physics such as virialization in the central region. 
Three-dimensional simulations (Kang et al. 1994; Navarro, 
Frenk, & White 1995; Evrard, Metzler, & Navarro 1995) 
showed that the gas is shock-heated to the virial tempera- 
ture and then settles into hydrostatic equilibrium (HSE) and 
that the shock separates the hydrostatic central region from 
the infalling flow. However, according to Navarro, Frenk, & 
White (1995), where iV-body/SPH simulations were used 
to examine the properties of x-ray clusters in the Q = 1 
CDM universe, the one-dimensional, self-similar solution of 
Bertschinger (1985) matches well the density and tempera- 
ture profiles of their simulated clusters except for the inner 
central region. Thus, the properties of accretion flow can be 
efficiently studied by one-dimensional calculations. 

The plan of the paper is as follows. In §2 we describe 
the numerical simulations. In §3 we present the results fo- 
cusing on the velocity of the accretion flows. In §4 we discuss 
the implications of the results on cosmology and on cosmic- 
ray acceleration by the accretion shocks. Details such as the 
evolution equations and the numerical scheme are given in 
the appendices. 

2 SIMULATIONS 

We follow the evolution of collisionless dark matter par- 
ticles accreting onto an initial density enhancement in an 
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otherwise homogeneous universe through calculations with 
a spherical iV-body code. The equations describing the par- 
ticle motion and gravity in comoving coordinates and the 
expansion of the background universe can be found, for ex- 
ample, in Peebles (1993) and are summarized in appendix A. 
The details of our numerical code are described in appendix 
B. 

The calculations start at the time corresponding to the 
initial expansion parameter = 10 -5 a o with 10 4 particles. 
Initially the particles are at rest and distributed uniformly, 
except for an excess mass at the origin to initiate the col- 
lapse. So the initial perturbation can be characterized by 
the point-mass perturbation 

5M(r) _(M(r)\ 1 

M(r) \ M J { ' 

where M(r) is the mass inside a radius r and M is a ref- 
erence mass. Note that the point-mass perturbation can 
be approximated by the constant power spectrum (i.e. , 
P(k) « costant with n = 0) (Hoffman & Shaham 1985). 
Although a scale- free power spectrum with n = — 1 would 
represent better the CDM power spectrum on the cluster 
scale (see e.g. , Bardeen et al. 1986 for the CDM power 
spectrum), the constant power spectrum would be a bet- 
ter approximation on scales a bit larger than that. Also ac- 
cording to Navarro, Frenk, & White (1995), the self-similar 
solution of Bertschinger (1985) with the above initial per- 
turbation represents well the structures of simulated clusters 
in an Q — 1 universe with the CDM power spectrum except 
for the inner central region. The amount of the excess mass, 
M , has been chosen so that at the present epoch about 1/3 
of the particles are placed inside the first caustic. It requires 
an excess mass corresponding to 1 — 2 particles. 

Note that the present value of the Hubble parameter, or 
h = H o /100 km s _1 Mpc~ 1 , enters the problem only through 
the normalization parameters such as t oc 1/h, L oc 1/h, 
and p oc h 2 (see the appendices). So, we can include the 
dependence on h implicitly by expressing the results in term 
of t h, L h, and ph~ 2 (or equivalently M h). 

If we had followed the evolution of BM with 7 = 5/3 
as well, its accretion shock would have occurred very close 
to the first caustic of DM and its accretion velocity out- 
side the shock would have been very similar to that of DM 
outside the first caustic. Since the accretion shock approxi- 
mately separates the inner virialized region from the outer 
infalling flow, here we will consider a cluster to be the region 
inside the first caustic. So, below, we define the radius of the 
first caustic, R c i, as the cluster radius, the mass inside R c i, 
M c i, as the total cluster mass, and the particle velocity just 
outside Rd, v acc , as the cluster accretion velocity. 

Our numerical set-up has been designed in a way paral- 
lel to that of Bertschinger (1985), so our numerical solution 
for the ft — 1 case should be identical to his self-similar solu- 
tion. Thus, the accuracy of our code has been tested against 
his solution. Figure 1 shows the present phase-space distri- 
bution of particles in the Einstein-de Sitter universe, which 
can be directly compared with Figure 6 of Bertschinger 
(1985). The variables used in the plot are related to those 
in the appendices as follows 



where rt a is the present turnaround radius. Our smooth- 
ing length in the gravitational force corresponds to A sm = 
0.0936 (see appendix B). The plot shows that the numerical 
solution for A J; A sm agrees well with the exact analytic so- 
lution. The position of the first caustic agrees well with that 
of the analytic solution and the accretion velocity outside it 
agrees almost exactly. On the other hand, for A <^ A sm the 
agreement in the two solutions is not so good. The numerical 
solution has given a smaller velocity than the analytic one 
and the positions of caustics have not been calculated cor- 
rectly, as expected. However, the quantities we are interested 
in, R c i, M c i, and v aC c, are less affected by the smoothing. 
Comparison of Figure 1 with Figure 6 of Bertschinger (1985) 
shows that these quantities have been calculated accurately 
within an error typically less than 5% or so. 

Figure 2 shows the evolution of the mass, radius, and 
accretion velocity of a cluster with the present mass of 
M ct (z = 0) = 10 15 /i _1 M Q as a function of redshift in 
the Einstein-de Sitter universe, in the open universe with 
Q = 0.2, and in the flat universe with Q = 0.2. The dotted 
lines represent the power-law evolution of the analytic solu- 
tion in the Einstein-de Sitter universe (Bertschinger 1985); 
Md oc (l+z)-\ Rd oc (l+z)- 4/:i , and v acc oc (l+z) 1/6 . The 
numerical solution in the Einstein-de Sitter universe follows 
the self-similar evolution very closely. Again, good agree- 
ment of the numerical solution with the analytic solution 
indicates that our code has been able to calculate the quan- 
tities used in this paper reliably. But those in other model 
universes do not show the power-law evolution in (1 + z), 
indicating the absence of the self-similarity. 



3 RESULTS 

Among the quantities we can extract from the x-ray obser- 
vations of clusters, the ICM temperature T x (usually aver- 
aged over a core region) is the one that can be estimated 
most reliably in hydrodynamic numerical simulations and 
that is less prone to numerical effects such as the problem of 
under-resolved structure. That is because T x is a conserved 
quantity per unit mass (i.e. the specific thermal energy). 
Although clusters have formed only recently (z < 1) and 
the dynamical time of cluster evolution is about the Hub- 
ble time, the structure of more-or-less relaxed clusters both 
in observations and simulations is well represented by the 
virialized region in HSE with nearly constant temperature 
(Forman & Jones 1982; Evrard, Metzler, & Navarro 1995). 
Exceptions are the clusters that have recently undergone 
mergers. Under the assumption of an isothermal sphere in 
HSE, the average temperature is proportional to M(< r)/r, 
or to the galaxy kinetic energy per unit mass, a 2 . So we 
use the ratio of the total cluster mass to radius, M c i/R c i, 
as a quantity to represent the observed, average temper- 
ature (emission weighted). As shown by three-dimensional 
simulations (Navarro, Frenk, White 1995; Evrard, Metzler, 
Navarro 1995), simulated clusters of different masses have 
similar structures when they are scaled with a radius of a 
fixed density contrast (S — p/p cr it ii 200). So the tempera- 
ture follows the scaling law, T oc M(< r c )/r c oc r 2 , where 
r c is the characteristic radius at a fixed density contrast. 
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Figure 1. Present phase-space distribution of dark matter particles from the numerical calculation in the Einstein-de Sitter universe. 
It has been drawn for direct comparison with Figure 6 of Bertschinger (1985). The relations of the variables used in the plot to those in 
the appendices are A = r/r ta and d\/d^ = (u + a r)(t /rt a ) — (8/9)A. Here, rt a is the present turnaround radius. 



Of course, the constant of proportionality for this scaling 
relation varies for different cosmologies. 

Here, we use the quantities, M c i, R c i, or M c i/R c i, to 
denote clusters. Note that R c i corresponds to the character- 
istic radius of S ~ 80. Under the set-up considered in this 
paper, only one of them is independent in a given cosmo- 
logical model universe, because of the similarity in cluster 
structure. According to Bertschinger (1985), the self-similar 
solution at the present epoch in the Einstein-de Sitter uni- 
verse gives 

^-2.22,-M P c( I5T ^ i -) 1/3 (4) 

T c! (0.3^) = 6.06 toV ( T5n ^) (6) 

v acc = 1.31 xIOVs^^) (7) 

Here R s is the shock radius, and T c i is the temperature at 
r = 0.3ii s , which we may take as an average ICM tem- 
perature of the self-similar flow. Navarro, Frenk, & White 
(1995) showed that the temperature profiles of simulated 
clusters normalized with the temperature at V2oo of 8 = 200 
can be represented within a factor of two by that of the 



self-similar solution of Bertschinger (1985) in r/r2oo ^ 0.3. 
In the core region within r/r2oo < 0.3, the temperature is 
approximately isothermal. Thus, T c i should be a reasonable 
approximation to the average temperature of the cluster core 
region. 

In Figure 3 we have plotted the present M c i/R c i for 
clusters with the present mass M ci {z = 0) = IO^/i-'Mq 
(upper panel) and for clusters with the present radius 
Rci(z = 0) = 2.5/i _1 Mpc (lower panel) in the open mod- 
els and flat models with fi D < 1 against the values of £l . 
Those with Q = 1 correspond to the values in the Einstein- 
de Sitter universe. Here we have chosen R c i = 2.5h~ Mpc as 
a fiducial value of the first caustic, since the typical radius 
of cluster core regions is often defined as r core — 0.5/i _1 Mpc 
and the Abell radius is Ra = 1.5/i _1 Mpc. According to the 
self-similar solution of Bertschinger (1985) given above, in 
the Einstein-de Sitter universe, 

Mcl i-.i . K, : : X1 • ( Md V /3 (8) 
R cl Mpc \l0^h-m o ) (8) 

for clusters with the given mass, while 

R cl - bJ6xW Mpc V2.5ft-iMpcJ (9) 

for clusters with the given radius. Comparison between these 
values and our numerical results for Q = 1 case shows about 
~ 10% error. Since M ct /R ct scales as M cl /R cl oc (M ct h) 2/3 
for clusters with different masses and M c i/R c i oc (R c ih) 2 
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Figure 2. Evolution of the radius, mass, and accretion velocity 
of a cluster with the present mass M c i(z = 0) = 10 15 h~ 'M0 
as a function of redshift in the Einstein-de Sitter universe (open 
squares), in the open universe with Q = 0.2 (filled circles), and 
in the flat universe with Q a = 0.2 (open circles). The dotted 
lines represent the power-law evolution in the exact self-similar 
solution in the Einstein-de Sitter universe (Bertschinger 1985). 



for clusters with different radii, the results can be scaled 
accordingly. 

The figure shows that the clusters, if they have the same 
mass or the same radius, are less tightly bound, and so have 
smaller M c i/R c i in the universe with smaller fl . Also, the 
clusters are less tightly bound in the flat universe with non- 
zero A than in the open universe, if £l is the same. This 
is due to the fact the background expands faster in the flat 
universe with non-zero A than in the open universe. In the 
open universe with Q a = 0.1, M c i/R c i is ~ 32% and ~ 68% 
smaller than that in the Einstein-de Sitter universe for clus- 
ters with given mass and radius, respectively. In the flat uni- 
verse with Q. = 0.1, M c i/Rci is ~ 40% and ~ 79% smaller. 
This means that clusters are less massive and less bound so 
the ICM temperature is lower in low density universes than 
in high density universes, if they are selected by a constant 
radius criterion. 

Figure 4 shows the value of M c i/R c i as a function 
of redshift in the Einstein-de Sitter universe, in the open 
universe with Q — 0.2, and in the flat universe with 
flo = 0.2. The upper panel is for clusters whose mass is 
M cl (z) = 1O 15 /i _1 M at given z. Similarly, the lower panel 
is for clusters whose radius is R c i(z) = 2.5/i _1 Mpc at given 
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Figure 3. Present M c i/R c i for clusters with the present mass 
M c i(z = 0) = 10 15 h~ 1 Mq (upper panel) and for cluster with 
the present radius R c i(z = 0) = 2.5h~ x Mpc (lower panel) in the 
open universes (filled circles) and the flat universes (open circles) 
with different Q a . That with f2 D = 1 represents the values in the 
Einstein-de Sitter universe. 



z. Thus, each point in the plots represent different clusters 
that should become heavier than 10 15 /i _1 Mq or bigger in ra- 
dius than 2.5/i _1 Mpc at the present epoch. In other words, 
the redshift dependences in these plots do not represent the 
evolution of a particular cluster in a "Lagrangian" sense. 
They tell us that clusters selected by either a constant mass 
or a constant radius criterion would be more tightly bound 
and so hotter in the past (at higher redshifts) than at the 
present epoch in all cosmologies. The dotted lines represent 
the power- law redshift dependences {M c i/R c i oc (1 + z) for 
clusters with a constant mass, and M c i/Rcl oc (1 + z) for 
clusters with a constant radius) of the self-similar solution 
in the Einstein-de Sitter universe. The numerical solution in 
the Einstein-de Sitter universe follows the self-similar solu- 
tions very well. Once again we can see that clusters are less 
tightly bound in lower density universes. The increases of 
Md/Rci with increasing redshift in both lower density uni- 
verses are slightly less than that in the Einstein-de Sitter 
universe. 

In Figure 5 we have plotted the accretion velocity, 
Vacc, for clusters with the present mass M c i(z = 0) = 
1O 15 /i _1 M0 (the first panel), for clusters with the present 
radius R c i(z = 0) = 2.5/i _1 Mpc (the second panel), and for 
clusters with the present M c i/R c i(z = 0) = 4x 10 14 M Q /Mpc 
(the third panel) in the open and flat universes with different 
o 's. Note that M ct /R a i = 4 x 1O 14 M /Mpc corresponds 
to T c i — 5.71 keV in the Einstein-de Sitter universe (see 
Eq. 6) (corresponding T c ; should be slightly lower in the 
low density universes). Clusters with different masses, radii, 
and Md/Ra's have the accretion velocity which scales as 
Vacc oc (M c ih) 1/3 , Vacc oc R c ih, and v aC c oc (M c i/R c i) 1/2 . 
Smaller accretion velocity in the universes with smaller 0, o 
is consistent with the fact that the clusters are less tightly 
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Figure 4. M c i/R c i as a function of redshift in the Einstein-de 
Sitter universe (open squares), in the open universe with f! = 0.2 
(filled circles), and in the flat universe with Q a = 0.2 (open cir- 
cles). The upper panel is for clusters with a mass at given z, 
M c i(z) = lO 15 h _1 M0. The lower panel is for clusters with a ra- 
dius at given z, R c i{z) = 2.5/i — 1 Mpc. The dotted lines represent 
the power-law evolution in the exact self-similar solution in the 
Einstein-de Sitter universe (Bertschinger 1985). 
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bound. The accretion velocity is slight smaller in the flat 
universe than in the open universe, if Q is same. The plot 
shows that the accretion velocity of clusters with a given 
mass or a given radius is smaller by up to ~ 41% and ~ 65% 
respectively in the flat universe with Q =0.1 than in the 
Einstein-de Sitter universe. 

For clusters with the same M c i / R c i , however, the differ- 
ence between different cosmological models is rather small. 
In the open universe with fi„ = 0.1 the accretion velocity 
is ~ 15% smaller than that in the Einstein-de Sitter uni- 
verse, while in the flat universe with Q =0.1 it is smaller 
by ~ 24%. This means we can expect that the clusters with 
the same observed temperature can have accretion velocities 
smaller only up to ~ 25% in low density universe depending 
on Q . 

Figure 6 shows the accretion velocity for clusters with 
a mass of M c i(z) — 1O 15 /i _1 M0 at given z (first panel), 
for clusters with a radius of R c i(z) — 2.5/i -1 Mpc at given 
z (second panel), and for clusters with M c i/R c i(z) = 4 x 
10 14 M Q Mpc at given z (third panel) in the Einstein-de Sit- 
ter universe, in the open universe with Q = 0.2, and in 
the flat universe with Q = 0.2. As in Fig. 4, they do not 
represent the evolutionary path of a cluster, but they show 
the dependence of v aC c on the redshift in a sample of clus- 
ters selected by a constant mass, or a constant radius, or a 
constant temperature criterion. The dotted lines represent 
the power-law dependence on the redshift of the self-similar 
solution in the Einstein-de Sitter universe, v a cc oc (1 + z) 1//2 
for clusters with a constant mass, v acc oc (1 + z) 3 ^ 2 for clus- 
ters with a constant radius, v aC c oc constant for clusters with 
a constant M c i/Rci- 



Figure 5. Present accretion velocity, v aC c, for clusters with the 
present mass M c ;(z = 0) = 10 15 /i — 'Mq (the first panel), for 
clusters with the present radius R c i(z = 0) = 2.5h~ x Mpc (the 
second panel), and for clusters with the present M c i/R c i = 4 X 
1O 14 M0/Mpc (third panel) in the open universes (filled circles) 
and the fiat universes (open circles) with different f! . That with 
Qo = 1 represents the values in the Einstein-de Sitter universe. 

The first two panels show that the accretion velocity of 
clusters with a constant mass or a constant radius was sig- 
nificantly larger in the past than now. On the other hand, 
the third panel shows that the accretion velocity of clus- 
ters with a constant M c i / R c i remains constant with redshift 
in the Einstein-de Sitter universe due to the scale-free na- 
ture. But it was larger at high redshifts in the model uni- 
verses with Q < 1. For instance, in the open universe with 
Q — 0.2 the accretion velocity was larger by ~ 8% at high 
redshifts, while in the flat universe with Q = 0.2 it was 
larger by ~ 18%. At high redshifts (z > 5) the difference 
in Vacc among different cosmological models considered here 
reduces to only a few percents. 



4 SUMMARY AND DISCUSSION 

In this paper, we have studied the properties of spheri- 
cally accreting flows onto an initially overdense perturba- 
tion in different cosmological model universes including the 
Einstein-de Sitter universe with Q = 1, the low-density 
open universe with Q < 1, and the low-density, flat uni- 
verse with Q < 1 and Qa = 1 — Qo- According to the 
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Figure 6. Accretion velocity, function of redshift in 

the Einstein-de Sitter universe (open squares), in the open uni- 
verse with Q a = 0.2 (filled circles), and in the flat universe with 
f2o = 0.2 (open circles). The first panel is for clusters with a 
mass at given z, M r j (2) = 10 15 h~ 'Mg. The second panel is 
for clusters with a radius at given z, R c i(z) = 2.5/i — 1 Mpc. The 
third panel is for clusters with Af c ;/i? c ; at given 2, M c i/R c i(z) = 
4 X 10 14 MqMpc. The dotted lines represent the power-law evo- 
lution in the exact self-similar solution in the Einstein-de Sitter 
universe (Bertschinger 1985). 



semi-analytic treatment of the self-similar solution for the 
Einstein-de Sitter universe by Bertschinger (1985), the ac- 
cretion shock of BM with 7 = 5/3 forms very close to the 
first caustic of DM, and the flow upstream to the accretion 
shock and the first caustic is identical regardless of whether 
it is collisional or collisionless. Thus, we have assumed that 
the position of the accretion shock can be approximated by 
that of the first caustic and the accretion velocity of BM out- 
side the accretion shock by that of DM outside the first caus- 
tic. Since self-similar solutions do not exist for Q 7^ 1 uni- 
verses, we have used a one-dimensional, spherical, iV-body 
code to study the properties of accreting matter. 

The accretion velocity onto a cluster with a given initial 
perturbation has decreased as the universe expands in all 
cosmological models. In the Einstein-de Sitter universe it 
has followed v aC c oc (1 + z) 1//6 (Bertschinger 1985), while in 
low density universes it decreases faster with time, especially 
at low redshifts (see Figure 2). This is related to the smaller 
deceleration of the expanding background. 

The properties of accretion flows around the clusters of 



the same characteristic (mass, radius, or temperature) will 
depend on the cosmological parameters, f2 and A, since 
the accretion rate is affected by the expansion rate of the 
background universe. The accretion velocity is smaller in 
a lower density universe and even smaller in a non-zero A 
universe. For clusters of a given mass at given redshift, the 
present accretion velocity in the low density universe with 
fi > 0.1 is smaller by up to ~ 41% than that in the Einstein- 
de Sitter universe. For clusters of a given radius at given 
redshift, the present accretion velocity in the low density 
universe with Q > 0.1 is smaller by up to ~ 65% than that 
in the Einstein-de Sitter universe. Hence, if the accretion 
velocity of infalling matter around clusters can be measured 
along with their mass and radius, then it could be used to 
discriminate the different cosmological models. 

However, for clusters of a given temperature or given 
Md/Rci, the present value of the accretion velocity as well 
as its evolution depend on the cosmological model rather 
weakly. According to the self-similar solution of Bertschinger 
(1985), in the Einstein-de Sitter universe, for clusters of 
a given temperature, the accretion velocity at present is 
given as v acc = 1.31 x 10 3 km s~ 1 (T ci /6.06keV) 1/2 , and it 
has been constant through the evolution of the universe. In 
the universes with smaller Q , the present accretion velocity 
is smaller by up to ~ 24% in the models considered here 
with 0.1 < f2 < 1. But this difference decreases to a few 
percents at high redshifts (see Figs. 5 and 6). Thus, it would 
be difficult to use the accretion velocity of clusters with a 
given temperature to discriminate among the cosmological 
models. 

However, this has an important implication for the 
model of the origin of ultra high energy cosmic rays that 
the acceleration of protons in accretion shocks around clus- 
ters can contribute significantly to the observed particle flux 
above 3 x 10 18 eV (Kang, Ryu, & Jones 1996; Kang, Rachen 
& Biermann 1996). We have shown in this paper that the 
matter accretion and the accretion shock around clusters are 
universal for all cosmological models, although the accretion 
regime is less strong in lower density universes. The shock 
velocity v a — (4/3)v aC c for the rich clusters of T c i — 10 keV, 
for example, is estimated to be 2200km s _1 in the Einstein- 
de Sitter universe, which is enough for the model to work. 
The results in this paper imply that the model might be ex- 
tended to all cosmologically viable universes of low density, 
since the reduction in the shock velocity is no more than 
~ 24%. 

The region of clusters observed by either optical or x-ray 
observations (r <J 1 — 1.5/i _1 Mpc) is inside the first caustic 
or the accretion shock {R c i ~ 1 — 3/i _1 Mpc for T c i ~ 1 — 10 
keV). Although the existence of hot gas heated by the ac- 
cretion shocks has been shown clearly in most cosmological 
simulations including hydrodynamics for any variants of cos- 
mological models (Kang et al. 1994; Cen & Ostriker 1994; 
Navarro, Frenk, & White 1995), any direct observations of 
gas near the shock would not be possible at present due to 
low surface brightness. But it has been suggested that the 
unbound hot gas of 10 5 — 10 6 K around clusters and groups 
heated by the large scale accretion shocks is a major com- 
ponent of the IGM which emits mostly at soft x-ray below 1 
keV (Ostriker & Cen 1996). This prediction might be tested 
by looking for the spatial correlation between the soft x-ray 
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cosmic background radiation and the observed large scale 
structure. 

On the other hand, the infalling clouds of unshocked, 
warm gas may be identified through the absorption lines of 
quasars which are located inside clusters of galaxies. This 
warm low density gas of 10 4 — 10 5 K is photoionized by the 
diffuse radiation from the hot postshock gas and the dif- 
fuse cosmic background radiation. In some studies (e.g. , 
Weymann et al. 1979; Foltz el al. 1986), the C IV absorp- 
tion systems of quasar emission lines with \z a bs — zqso\ < 
3000km s _1 are interpreted as clouds associated with rich 
clusters where the quasars reside. The characteristics of 
these systems of C IV absorbers are different from those of 
the typical intervening C IV absorbers. It was noted that the 
velocity difference is somewhat too large compared to the 
typical velocity dispersion of galaxies (400 — 1200 kms -1 ) 
in rich clusters (Foltz el al. 1986). The accretion velocity, 
however, is a bit larger than the galaxy velocity dispersions, 
since it is given by v acc = 1.31 x 10 3 km s _1 (T ci /6.06keV) 1/2 
in the Einstein-de Sitter universe. Thus it is possible that 
these absorption systems are in fact the infalling clumps of 
gas upstream of the accretion shock. 

In follow-up papers we will study the properties and 
statistics of three-dimensional accretion flows in simulated 
universes by analyzing the data from three-dimensional hy- 
drodynamic simulations of various cosmological models (e.g. 
, Kang el al. 1994; Cen & Ostriker 1994) . These simulations 
have been performed by a high-resolution, grid-based, Eu- 
lerian code (Ryu et al. 1993) which resolves the shock dis- 
continuity in 2-3 cells and is designed specifically to handle 
the flows with supersonic bulk motions. Thus, the low den- 
sity regions around caustics and shocks are well represented 
in these simulations, even though the high density core re- 
gion of clusters might be under-resolved (Kang et al. 1994). 
The present one-dimensional study will provide some useful 
guidance for such studies. 
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APPENDIX A: BASIC EQUATIONS 

In comoving coordinates, the equations to describe the evo- 
lution of the collisionless dark matter particles are written 

as 



^4nG[p(r)-p], 



(A3) 



dr 1 

777 = _u ' 

at a 

and 

du a 1 _ , 

— = u- -V(j>, 

at a a 

where r is the comoving position, u is the proper peculiar 
velocity, cj> is the proper peculiar gravitational potential, and 
a is the cosmic expansion parameter. The potential is given 
by the Poisson equation 



(Al) 



(A2) 



where p(r) is the comoving matter density, p is the average 
comoving matter density, and a is the present value of the 
the expansion parameter. 

The expansion parameter and its time derivative (or the 
expansion rate) are calculated as a function of time by the 
following equations (see Peebles 1993). 



a 
a 
a 
a„ 



2/3 



■ (H 

2 (1 - 
H t = 



for Einstein — deSitter (A4) 



(cosh?) — 1) and 



2(1 -Q 

1/3 



.3/2 



(sinh ij — rj) , 



for open (A5) 



i = (l^fc) Sinh2/3 {l^^H ot ) ^ for flat (A6) 



and 



(A7) 



Here, H is the present value of the Hubble parameter and 
Q. , Qr, and Qa are constants. The density parameter is 
given as 

StvGp 



Q. 



(A8) 



which is the present average mass density in terms of the 
critical density. The parameter associated with the radius 
of curvature R is given as 

9 (A9) 

(a H R) 



which is positive for an open universe and zero for others. 
The parameter associated with the cosmological parameter 
A is given as 

n A = ^ (Aio) 

which is positive for a flat universe and zero for others. 
The three parameters give the relative contributions to the 
present expansion rate by satisfying 



n + n R + n A 



(All) 



We have solved numerically the equations governing the 
evolution of the dark matter particles [Eqs. (A1)-(A3)], si- 
multaneously with the equations describing the background 
universe [Eqs. (A4)-(A7)]. In the next appendix, we describe 
the scheme used. 



APPENDIX B: NUMERICAL SCHEME 

In the one-dimensional spherical code to calculate the time 
integration of the equations of motion [Eqs. (A1)-(A2)], 
we have adopted the second-order accurate Lax-Wendroff 
scheme instead of the more popular leapfrog scheme, since 
the code was originally designed to be a part of a one- 
dimensional code which follows the evolution of the baryonic 
matter as well as that of the dark matter in a way parallel 
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to the three-dimensional cosmological hydrodynamic code 
(Ryu et al 1993). So the update of the position and veloc- 
ity of t he particle i from n time step to n + 1 is carried out 



by the following two steps: 

n + 1 /2 _ n , Mi 

r ' ~ Tl + 2 a" ' 



1+1/2 



„ At"d' 1 „ At n B 



(Bl) 
(B2) 



and 



n+1 



Tl I Ail 

= n +At 



IT 



i±U2 



At n 



a n + l/2 ' 
,n+l/2 



(i 



-1/2 



'+1/2 



At 



n n + 1/2 



where r denotes the gravitational force. With the variables 



normalized with the present age of the universe, t ci the co- 



moving l size of the box. Z/n. and o, ql 1 is calculated with 



-7T~ (H t ) 2 Q 
8tt (a n ) 2 



al SM(rf) 



OP 



(B5) 



1/2 



where <5M(r") is the mass excess inside r™. Similarly, p. 
is calculated with the quantities at n + 1/2. Here e is the 
smoothing parameter which prevents the time step becom- 
ing too short, or prevents the particles being accelerated 
anomalously for a given time step around the origin. 

A particle, which is placed close to and approaches the 
origin, can pass the origin in the half time step (n + 1/2). 
Then, the absolute value of its position is used to calculate 
the gravitational force at the half time step. If a particle 
passes the origin after the full time step (n+1), its position 
and velocity are reset to their negative values. 

The time step, At, is determined so that Aa/a is con- 
stant in each time step. In all the calculations discussed in 
this paper, we have used 10 4 particles and Aa/a — 10 -3 . 
With these, e = 5 x 10~ 2 in unit of L assures that the grav- 
itational time scale in the core region with r < e is com- 
fortably small compared to the time step along the whole 
calculation. However, the force smoothing makes the den- 
sity smoothed in the core region with r <^ e (see §2). 
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